function [h6]=channel_plot(h)
  
 rate=1;
 h2=h(:,1:3:end);
 %升采样

% 复制法
    h3=h2(1:16,:);
    h4=repmat(h3(1,:),512,1);
    for i = 2:size(h3,1)
        h4=[h4;repmat(h3(i,:),512,1)];
    end
    h5=h4/max(max(abs(h4)));
    h6=h5;
    h6(abs(h6)<rate*mean(mean(abs(h6))))=0;
%  插值法
    h3=h2(1:17,:);
    % xx=1:1/512:17-1/512;
    % h4=interp1(1:17,h3,xx,'spline');
    % h5=h4/max(max(abs(h4)));
    % h6=h5;
    % h6(abs(h6)<rate*mean(mean(abs(h6))))=0;

% 什么也不做
% h6=h/max(max(abs(h)));

[M,N]=size(h6);
Fn = dftmtx(N);       % Generate the DFT matrix
Fn = Fn./norm(Fn);    % normalize the DFT matrix

Fm = dftmtx(M);       % Generate the DFT matrix
Fm = Fm./norm(Fm);    % normalize the DFT matrix   
    %画信道图
% 时间-时延
% 时间
% |
% |
% |-------->时延
    fb_tao=3/16000;  %s
    fb_t=510*4/16000/512;  %s

    t=(0:(M-1))*fb_t;  
    doppler=-1/fb_t/2:1/fb_t/(M-1):1/fb_t/2;
    tao=(0:(N-1))*fb_tao;
    freqence=-1/fb_tao/2:1/fb_tao/(N-1):1/fb_tao/2;
    hlog=20*log10(abs(h6));
    set(figure,'visible','on');
    figure(1)
    pcolor(tao*1000,t,hlog);
%     mesh(0:fb_tao:85*fb_tao,0:fb_t:8191*fb_t,hlog);
    colormap(turbo);
    shading interp;
    tao_min=min(tao*1000);tao_max=max(tao*1000);%ms
    t_min=min(t);t_max=max(t);
    c_min=-40; c_max=0;clim([c_min c_max]);
    axis([tao_min,tao_max,t_min,t_max]);   
    xlabel( 'delay (ms)','FontName','Times New Roman','FontSize',12 );
    ylabel( 'Time(s)','FontName','Times New Roman','FontSize',12 );
    title('Time evolution (dB)','FontName','Times New Roman','FontSize',12 );
% 时间-频率
% 时间
% |
% |
% |-------->频率
    h7=fftshift(h6*Fn,2);
    hlog=20*log10(abs(h7));
    set(figure,'visible','on');
    figure(2);
%   pcolor(0:1/fb_tao/85:1/fb_tao,0:fb_t:8191*fb_t,hlog);
    mesh(freqence,t,hlog);
    colormap(turbo);
    shading interp;
    xlabel( 'frequence (Hz)','FontName','Times New Roman','FontSize',12 );
    ylabel( 'Time(s)','FontName','Times New Roman','FontSize',12 );
    title('Time evolution (dB)','FontName','Times New Roman','FontSize',12 );

    % 多普勒-时延
% 多普勒
% |
% |
% |-------->时延
 
    h8=fftshift(Fm*h6,1);
    h9=h8.*conj(h8);
    hlog=10*log10(h9./max(max(h9)));

    set(figure,'visible','on');
    figure(3)
    pcolor(tao*1000,doppler,hlog);
%   mesh(0:fb_tao:85*fb_tao,0:1/fb_t/8191:1/fb_t,hlog);;
    colormap(turbo);
    shading interp;
    c_min=-40; c_max=0;clim([c_min c_max]);%colormap( hot );
    tao_min=min(tao*1000);tao_max=max(tao*1000);t_min=min(doppler);t_max=max(doppler);
    axis([tao_min,tao_max,t_min,t_max]);
    xlabel( 'delay (ms)','FontName','Times New Roman','FontSize',12 );
    ylabel( 'Doppler(Hz)','FontName','Times New Roman','FontSize',12 );
    title('Time evolution (dB)','FontName','Times New Roman','FontSize',12 );
    